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EAE FORTRAN 


These few paragraphs summarize the revisions made to the PDP-8 
FORTRAN operating system to utilize the extended arithmetic unit option 
(Type 182 EAE). Four arithmetic routines were rewritten—alignment, 
normalize, multiply and divide. The reduction in execution time that was 
achieved was quite significant. 

For example, an iterative-type problem that used several function 
subroutines (sin, cos, etc, ) had been averaging 35 seconds per printout 
with non-EAE FORTRAN, This execution time was cut to 17 seconds with 
the EAE FORTRAN package. Comparative execution times of. the subroutines 
for non-EAE FORTRAN versus EAE FORTRAN are shown in Table I. 

Another improvement besides the faster execution time was gained 
with EAE FORTRAN. Since the multiply routine calculates a full 48-bit 
product and rounds instead of truncates to 24 bits, an increase in significance 
of the product was noted. For example, starting with 0,1 and successively 
multiplying by 10 for 100 times should yield 0.lOOOOOE+100. Non-EAE FORTRAN 
yielded 0,999973E+99, while EAE FORTRAN yielded 0.999999E+99, 

These modifications work with the FORTRAN operating system of 
March 2, 1967. They have not been tested with any other version, but would 
^'probably" work. No changes must be made in operating procedure or any 
other portion of the program, as this modification loads over the regular 
arithmetic subroutines. 


This word was performed under the auspices of the U. S. Atomic Energy Commission. 
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TABLE I 

Comparative Execution Times of EAE and Non-EAE Arithmetic Subroutines 


Routine Execution Time (usee) 



Minimum 


Average 

Maximum 


EAE 

Non-EAE 

EAE 

Non-EAE 

EAE 

Non-EAE 

Alignment 

119 

90 

125 

423 

130 

750 

Normalize 

41 

102 

51 

469 

61 

836 

Multiply 

239 

1116 

303 

1160 

366 

1204 

Divide 

128 

1893 

183 

1936 

237 

1978 
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/EAE FORTRAN MODIFICATIONS 




*103 



0]03 

2224 

ALGN 





*141 



0 ] 4 ] 

2322 

MP2 



0]42 

2320 

MP4 





/NORMALIZE ROUTINE FOR 

EAE FORTRAN 



*2200 



2200 

0000 

NORM, 

0 


2201 

7300 


CLA CLL 


2202 

1445 


TAD 1 L2 


2203 

7450 


SNA 

/IS LEAST SIG. 0? 

2204 

5220 


JMP ZCHECK 

/YES, CHECK MOST SIG. 

2205 

7421 

Bl, 

MOL 

/LEAST SIG. TO MO 

2206 

1444 


TAD 1 H2 

/MOST SIG. TO AC 

2207 

7411 


NMI 

/NORMALIZE 

2210 

3444 


DCA 1 H2 

/DEPOSIT INTO 2ND ARG. 

2211 

7501 


MQA 


2212 

3445 


DCA 1 L2 


2213 

7641 


SCA CLA 

/GET STEP COUNTER 

2214 

7041 


CIA 

/NEGATE 

2215 

1443 


TAD 1 X2 

/ADD EXPONENT 

2216 

3443 

B2, 

DCA 1 X2 


2217 

5600 


JMP 1 NORM 

/EXIT 

22?0 

1444 

ZCHECK 

TAD 1 H2 


2221 

7640 


SZA CLA 

/MOST SIG. 0? 

2222 

5205 


JMP Bl 


2223 

5216 


JMP B2 

/YES, 0 TO EXPONENT 



/ALIGNMENT ROUTINE FOR 

EAE FORTRAN 

2224 

0000 

ALGN, 

0 


2225 

7300 


CLA CLL 


2226 

3303 


DCA ID 


2227 

1440 


TAD 1 XI 


2230 

7041 


CIA 


2231 

1443 


TAD 1 X2 


2232 

7450 


SNA 

/EXPONENTS EQUAL? 

2233 

52-75 


JMP DONE 

/YES 

2234 

7510 


SPA 

/X2 GREATER 

2235 

5246 


JMP TWO 









2236 

7041 


CIA 

•2237 

7040 


CMA 

2240 

3271 


DCA SHFT 

2241 

7340 


CLA CLL CMA 

2242 

3303 


DCA ID 

2243 

1041 


TAD HI 

2244 

3304 


DCA PNT 

2245 

5254 


JMP BOTH 

2246 

7040 

TWO 

CMA 

2247 

3271 


DCA SHFT 

2250 

1044 


TAD H2 

2251 

3304 


DCA PNT 

2252 

1440 


TAD 1 XI 

2253 

3443 


DCA 1 X2 

2254 

1271 

BOTH, 

TAD SHFT 

2255 

1306 


TAD M24B 

2256 

7700 


SMA CLA 

2257 

5300 


JMP NOGO 

2260 

7101 


CLL lAC 

2261 

1304 


TAD PNT 

2262 

3305 


DCA PNT+1 

2263 

1705 


TAD 1 PNT+1 

2264 

7421 


MOL 

2265 

1704 


TAD 1 PNT 

2266 

7510 


SPA 

2267 

7020 


CML 

2270 

7415 


ASR 

2271 

0000 

SHFT, 

0 

2272 

3704 


DCA 1 PNT 

2273 

7501 


MQA 

2274 

3705 


DCA 1 PNT+1 

2275 

7300 

DONE, 

CLA CLL 

2276 

2224 


ISZ ALGN 

2277 

5624 


JMP 1 ALGN 

2300 

2303 

NOGO, 

ISZ ID 

2301 

4561 


JMS 1 SWP 

2302 

5624 


JMP 1 ALGN 

2303 

0000 

ID, 

0 

2304 

0000 

PNT, 

0 

2305 

0000 


0 

2306 

7750 

M24B, 

-30 


/X2-X1-1 

/NUMBER OF SHIFTS 


/SET POINTER 
/X1-X2-1 

/NUMBER OF SHIFTS 
/SET POINTER 
/X2=X1 


/TOO MANY SHIFTS? 
/YES 


/LEAST SIG. TO MQ 
/MOST SIG. TO AC 
/POSITIVE? 

/NO, SET MINUS SIGN 
/SHIFT RIGHT 


/EXIT 

/MAKE ARG2 LARGEST 
/EXIT 
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/DIVIDE ROUTINE FOR EAE FORTRAN 


/CALLING ROUTINE MODIFICATIONS 


205 ] 

7000 

*2051 

NOP 

2063 

1005 

*2063 

TAD SIGN 

2064 

7710 


SPA CLA 

2153 

1005 

*2153 

TAD SIGN 

2154 

7700 


SMA CLA 

2155 

5500 


JMP1 INT 1 

2156 

1444 


TAD 1 H2 

2157 

7041 


CIA 

216jel 

3444 


DCA 1 H2 

2161 

5500 


JMP1 INT 1 

2400 

0000 

*2400 

DDVD, 

0 

2401 

4672 


JMS 1 SN 

2402 

3006 


DCA DIVID 

2403 

1442 


TAD 1 LI 

2404 

7104 


CLL RAL 

2405 

3224 


DCA DDVD3 

2406 

1441 


TAD 1 HI 

2407 

7004 


RAL 

2410 

3245 


DCA DDVD4 

2411 

1245 


TAD DDVD4 

2412 

3217 


DCA DDVD2 

2413 

1445 


TAD 1 L2 

2414 

7421 


MOL 

2415 

1444 


TAD 1 H2 

2416 

7407 


DVl 

2417 

0000 

DDVD2 

0 

2420 

3217 


DCA DDVD2 

2421 

7501 


MQA 

2422 

3444 


DCA 1 H2 

2423 

7405 


MUY 

2424 

0000 

DDVD3 

0 

2425 

7141 


CLL CMA lAC 

2426 

1217 


TAD DDVD2 

2427 

7450 


SNA 

243^ 

5673 


JMP 1 DV5 1 

2431 

7420 


SNL 

2432 

5250 


JMP DDVD6 


5 


/CHECK SIGNS 
/CLEAR FLAG 


/SHIFT DIVISOR LEFT 


/LEAST SIG. TOMQ 
/MOST SIG. TO AC 
/(MS2+LS2*2**-12)/MS1 

/REMAINDER 

/QUOTIENT 

/QUO.*LS1*2**-12 
/NEGATE PRODUCT 

/HI PROD.=REM. (L=l) 
/HI PROD.>REM. 













2433 

3224 


DCA DDVD3 

/REM.>HI PROD. 

2434 

7040 


CMA 


2435 

3006 


DCA DIVIO 

/SET FLAG 

2436 

750 ] 


MQA 

/SUBTRACT LOW ORDER PROD. 

2437 

7141 


CLL CIA 

/FROM DDVD3+0*2**-12 

2440 

7421 


MQL 


2441 

7420 


SNL 


2442 

7040 


CMA 

/DDVD3-1 

2-443 

1224 


TAD DDVD3 


2444 

7407 


DVI 

/(REM.-QUO.*LSl*2**-12)/MSl 

2445 

0000 

DDVD4, 

0 


2-446 

5647 


JMP 1 .+1 


2-447 

27^1 


AAADIV 


2450 

704 ] 

DDVD6, 

CIA 


2451 

3217 


DCA DDVD2 

/HI PROD.>REM. 

2452 

1245 


TAD DDVD4 


2453 

7141 


CLL CIA 


2454 

1217 


TAD DDVD2 


2455 

7420 


SNL 


2456 

5263 


JMP .+5 


2-457 

3217 


DCA DDVD2 


2460 

7040 


CMA 


2461 

1444 


TAD 1 H2 


2462 

3444 


DCA 1 H2 


2463 

7200 


CLA 


2464 

1217 


TAD DDVD2 


2465 

744j0 


SZA 


2466 

5244 


JMP DDVD4-1 


2467 

567J0 


JMP 1 DV5M1 


2470 

27J04 

DV5M1, 

DDVD5-1 


2471 

56J010 

RDV, 

JMP1DDVD 


2472 

2664 

SN, 

SGN 


2473 

2705 

DV5I, 

DDVD5 




*2701 



27^1 

7701 

AAADIV, 

MQA CLA 

- 

2702 

7120 


CLL CML 

/FOR REM.>HI PROD. CASES 

2703 

2006 


ISZ DIVIO 

/REM>HI PROD.? 

2704 

7141 


CLL CIA 

/NO 

2705 

3445 

DDVD5, 

DCA 1 L2 


2706 

7420 


SNL 


2707 

7040 


CMA 

/H.O. QUO.-l 

27 ] 0 

1444 


TAD 1 H2 


2711 

7510 


SPA 


2712 

5716 


JMP 1 DV7I 


2713 

3444 


DCA 1 H2 


2714 

5715 


JMP 1 RET 


2715 

2471 

RET, 

RDV 


2716 

2307 

DV7I, 

DDVD7 









23j2f7 

7110 

*2307 

DDVD7, 

CLL RAR 

2310 

3444 


DCA 1 H2 

2311 

1445 


TAD 1 L2 

2312 

7010 


RAR 

2313 

3445 


DCA 1 L2 

2314 

2443 


ISZ 1 X2 

2315 

5717 


JMP 1 RETl 

2316 

5717 


JMP 1 RETl 

2317 

2471 

RETl, 

RDB 


/SHIFT QUO. RIGHT 
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/MULTIPLY ROUTINE FOR EAE FORTRAN 


*2600 


2600 

0000 

OMUL, 

0 


2601 

4264 


JMS SON 

/CHECK SIGNS 

2602 

1064 


TAD M2 


2603 

3361 


DCA CNTl 


2604 

1254 


TAD BSTKAD 


2605 

3017 


DCA TEMC 

/PRODUCT TABLE POINTER 

2606 

7040 


CMA 


2607 

1041 


TAD HI 


2610 

3015 


DCA TEMA 

/ARGl POINTER 

2611 

1064 

A2LO, 

TAD M2 


2612 

3362 


DCA CNT2 


2613 

7040 


CMA 


2614 

1044 


TAD H2 


2615 

3016 


DCA TEMB 

/ARG2 POINTER 

2616 

1415 


TAD 1 TEMA 

/GET ARGl 

2617 

3222 


DCA MUL 


2620 

1416 

A3LO, 

TAD 1 TEMB 

/GET ARG2 

2621 

7425 


MOL MUY 


2622 

0000 

MUL, 

0 


2623 

3417 


DCA 1 TEMC 

/STORE IN PRODUCT TABLE 

2624 

7501 


MQA 


2625 

3417 


DCA 1 TEMC 


2626 

2362 


ISZ CNT2 


2627 

5220 


JMP A3LO 


2630 

2361 


ISZ CNTl 


2631 

5211 


JMP A2LO 


2632 

1263 


TAD A 

/ADD PRODUCTS FROM TABLE 

2633 

1262 


TAD B 


2634 

1260 


TAD D 


2635 

7004 


RAL 

/CARRY 

2636 

7630 


SZL CLA 


2637 

7001 


lAC 

/ROUND 

2640 

7100 


CLL 


2641 

1261 


TAD C 


2642 

1257 


TAD E 


2643 

1256 


TAD F 


2644 

3445 


DCA 1 L2 


2645 

7004 


RAL 


2646 

1255 


TAD G 


2647 

3444 


DCA 1 H2 


2650 

1005 


TAD SIGN 


2651 

7710 


SPA CLA 


2652 

4537 


JMS 1 NEG 


2653 

5600 


JMP 1 DMUL 


















2654 

2654 

BSTKAD, 

G-1 

2655 

0000 

G, 

0 

2656 

0000 

F, 

0 

2657 

0000 

E, 

0 

266jZf 

0000 

D, 

0 

2661 

0000 

c. 

0 

2662 

0000 

B, 

0 

2663 

0000 

A, 

0 

2664 

0000 

SGN, 

0 

2665 

1444 


TAD 1 H2 

2666 

0300 


AND MASK 

2667 

1441 


TAD 1 HI 

2670 

3005 


DCA SIGN /SIGN OF PROD. OR QUO. 

2671 

1444 


TAD 1 H2 

2672 

771,0 


SPA CLA 

2673 

4537 


JMS 1 NEG 

2674 

1441 


TAD 1 HI 

2675 

77] 0 


SPA CLA 

2676 

4536 


JMS 1 NEGl 

2677 

5664 


JMP 1 SGN 

2700 

4000 

MASK, 

*2761 

4000 

2761 

0000 

CNTl, 

0 

2762 

0000 

CNT2, 0 

/MULTIPLY FOR XTRA-SPECIAL INDEX CALCULATOR 

*2320 

2320 

0000 

MP4, 

0 

2321 

7425 


MOL MUY 

2322 

0000 

MP2, 

0 

2323 

5720 


JMP 1 MP4 
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/ DEFINITIONS FOR EAE FORTRAN 


X\=40 

HI =41 

LI =42 

X2--43 

H2-44 

L2-45 

SIGN=5 

TEMA=15 

TEMB=16 

TEMC=17 

M2^-64 

SWP=161 

I NT 1=100 

DIVID=6 

NEG=137 

NEG1=136 

DV 1=7407 

NM 1=7411 

SHL-7413 

ASR=7415 

MQL=7421 

MUY=7405 

MQA=7501 

SCA=7441 














t 




A 

2663 

NEG 

0137 

AAADIV 

270] 

NEGl 

0136 

ALGN 

2224 

NMI 

7411 

ASR 

7415 

NOGO 

2300 

A2LO 

2610 

NORM 

2200 

A3LO 

2616 

PNT 

2304 

B 

2662 

RDV 

2471 

BOTH 

2254 

RET 

2715 

BSTKAD 

2654 

RETl 

2317 

B1 

2205 

SCA 

7441 

B2 

2216 

SGN 

2664 

C 

2661 

SHFT 

2271 

CNT1 

2652 

SHL 

7413 

CNT2 

2653 

SIGN 

0005 

D 

2660 

SN 

2472 

DDVD 

2400 

SWP 

0161 

DDVD2 

2417 

TEMA 

0015 

DDVD3 

2424 

TEMB 

0016 

DDVD4 

2445 

TEMC 

0017 

DDVD5 

2705 

TWO 

2246 

DDVD6 

2450 

XI 

0040 

DDVD7 

2307 

X2 

0043 

DIVIO 

0006 

ZCHECK 

2220 

DMUL 

2600 



DONE 

2275 



DVI 

7407 



DV5I 

2473 



DV5M1 

2470 



DV7I 

2716 



E 

2657 



F 

2656 



G 

2655 



HI 

0041 



H2 

0044 



ID 

2303 



INTI 

0100 



LI 

0042 



L2 

0045 



MASK 

2700 



MP2 

2322 



MP4 

2320 



MQA 

7501 



MOL 

7421 



MUL 

2620 



MUY 

7405 



M2 

0064 



M24B 

2306 




n 











